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We propose a quantum description of the cooling of a micromechanical flexural oscillator by a 
one-dimensional transmission line resonator via a force that resembles cavity radiation pressure. 
The mechanical oscillator is capacitively coupled to the central conductor of the transmission line 
resonator. At the optimal point, the micromechanical oscillator can be cooled close to the ground 
state, and the cooling can be measured by homodyne detection of the output microwave signal. 
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I. INTRODUCTION 



Micro- and nano-mechanical resonators have been an 
interesting research topic due to their broad application 
in technology and fundamental physics [1| . This includes 
studies of ultrahigh precision displacement detection ,2] , 
mass detection [3], gravitational- wave detectors [1, 
and attempts to observe quantum behavior of mechan- 
ical motion [i, 0, B i, Many of the applications 
are fundamentally limited by thermal fluctuations, and 
in order to reduce their effects, it is desirable to cool 
the mechanical oscillators. Recently, various schemes 
like the laser sideband cooling schemes developed for 
trapped ions and atoms llj, have been proposed for 
significantly cooling a mechanical resonator (MR) cou- 
pled to a Cooper-pair box [H, [H, [13, [Si ) a flux qubit 
[TgI [17| , a superconducting single-electron transistor [l^ , 
quantum do ts Il9l, trapped ions [201 , and optical cavi- 
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[211 122, I23UM [2^, |27|, |28|, |29|, |30|, |31|, M, M, 
3a . [39| . On the experimental side, op- 
tomechanical cooling schemes have been shown to be 
promising [H [H, [H, H, [H, [H, [H [H : the MR was 
cooled to ultra-low temperatures via either photothermal 
forces or radiation pressure by coupling it to a driven cav- 
ity. There are two main optomechanical cooling schemes. 
The first one involves an active feedback loop [2J, [2^, Is^l , 
and the second one works via passive back-action cool- 
ing (also called self-cooling) [22, 2l [H [2|. A fully 
quantum-mechanical description of cavity-assisted cool- 
ing schemes for optomechanical systems has been give n 
in Refs. [H, [H, [SJ, [H, [H, (for a review, see [38ll39t). 
Ground-state cooling of a mechanical resonator via pas- 
sive cooling schemes based on radiation pressure has also 
been investigated theoretically [H, [sl, HI, ■ 

Recently, other optomechanical-like cooling schemes 
have been proposed to replace the optical cavity by a 
radio-frequency (RF) circuit 40j,|4l| or a one-dimensional 
transmission line resonator (TLR) However, the 

theoretical understanding of the cooling schemes via a 



RF circuit in Refs. [40|, |4l| or via a TLR in Ref. [43 

is based on a classical description of the motion of the 
MR. A quantum-mechanical description of the motion 
of the MR, in a similar system consisting of a mechani- 
cal resonator capacitively coupled to a superconducting 
coplanar waveguide (which is an example of a TLR) , was 
discussed recently in Ref. [11] , which focused on studying 
the entanglement between the MR and the TLR without 
considering the cooling of MR. Most recently, Teufel et al. 
Q considered the cooling of a MR by applying directly 
the theoretical analysis of the cavity-assisted back-action 
cooling scheme [s^ to a superconducting microwave res- 
onator. They also presented experimental data about 
the cooling effect on the MR due to the microwave radia- 
tion field. The quantum-mechanical description of TLR- 
assisted cooling of a MR has also been investigated in 
Ref. [i^ via embedding a SQUID [46], which allows to 
control the coupling strength between MR and TLR by 
controlling the flux through the SQUID. 

There are some practical advantages [13, I3 
microwave TLR schemes. The TLR is realized in a thin 
on-chip superconducting film and is easily pre-cooled by 
standard dilution refrigeration techniques. It is ready to 
be integrated with quantum circuits containing Joseph- 
son junctions which may offer a sensitive measurement 
and a connection with quantum information processing. 
In addition, the size of the mechanical resonator could 
be much smaller than the wavelength of the radiation in 
the TLR, unlike in optical cavity experiments that work 
with refiection. 

In this paper, we present a quantum-mechanical de- 
scription and use it to investigate the motion of the 
MR when it is coupled capacitively to a driven TLR 
as in Ref. ^43| where the calculation was carried out 
in a semi-classical framework. The Hamiltonian of our 
TLR-assisted model is also studied in Refs. [i^ |4^ . 
and is very similar to that of a MR coupled to a 
driven optical cavit y v ia radiation pressure coupling 
P, m, m HE m, M m, HI. we study the TLR- 
assisted passive back-action cooling of a MR in detail by 
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using a quantum Langevin description (without taking 
into account quantum entanglement [43] between the MR 
and the TLR). One of the main results of our work is to 
show that the MR can be cooled close to its ground state 
using realistic parameters: final effective mean phonon 
numbers below 1 can be reached assuming an initial tem- 
perature of 10 mK which can be achieved using a dilu- 
tion refrigerator. We discuss in detail how such a ground 
state cooling of the MR can be obtained for all kind of 
parameter choices in practice. 



II. MODEL AND HAMILTONIAN 

The system that we consider is shown schematically in 
Fig.[TJ a MR is fixed on both ends (or a cantilever fixed on 
one end) located at the center of the TLR and is coupled 
capacitively to the central conductor of the TLR [49[ . We 
restrict the description to the fundamental flexural mode 
of oscillation of the MR which is modeled as a harmonic 
oscillator of frequency LOh and effective mass m. The TLR 
is driven by an external microwave at a frequency LOd 
and can be modeled as a single mode LC resonator with 
frequency = 1/ \/LaCa (the second mode of the TLR 
[soj). where La is the inductance and Ca the capacitance 
of the TLR. 
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FIG. 1; (Color online) Schematic diagram of a mechanical 
resonator (MR) located at the center of a one-dimensional 
transmission- line resonator (TLR). The external microwave 
drive field enters from the left and drives the TLR. The sig- 
nal at the output on the right end can be used to measure the 
motion of the MR via homodyne detection [13, H, 113, HJ] . 



the z-direction [d is the initial equilibrium distance with- 
out the coupling and Cg the corresponding initial capaci- 
tance; typically d ~ 1 /xm [13]). The capacitor is assumed 
to be placed in the center of the structure, i.e., its volt- 
age is given by the antinode voltage of the second mode: 
V = Vrmsia'' + a) (where Vrms = \J^'-^'al^a. is the rms 
voltage i49]), since the length of the MR is usually much 
shorter than that of the TLR: i cm > Z 10 - 100 
/im. The last term in Eq. ([1]) describes the input driving 
of the TLR by an external microwave field with the cou- 
pling strength [H [H [H |£| = ^^nPjhw'a, where k is 
the decay rate of the TLR, P is the input external mi- 
crowave drive power. Here, the non-rotating wave terms 
like ae"*'^''* and ate*""** have been ignored since we keep 

kl < < 

Usually, the fundamental oscillation frequency is of the 
order of 27rx (10"^ - 10^) Hz for micromechanical res- 
onators and 27rx (10^ - 10^) Hz for nanomechanical res- 
onators; the TLR frequency can be made to be of the 
order of 27r x 10 GHz. Here we will focus on the case of 
a micro-MR for which the condition ^ lo'^ is satisfied. 
In the interaction picture with respect to Hujda^a and 
neglecting the rapidly-oscillating terms, the Hamiltonian 
reads 



Hi = hAoa^a + [ — 
2m 



b ^2 



%0 



(2ata + l)x + h{ea^ + e*a) 



(2) 



where go ■— CgVj^^g/{hd) is a real coupling constant; 
Aq = oja — oJd is the detuning, and uja = uj'a + C^V^^g/h 
is the modified frequency of the TLR shifted by the cou- 
pling between TLR and MR. 

This Hamiltonian resembles that used in cavity- 
assisted cooling schemes [H, [H, [sj, [H, [H, [s^l- This 
suggests that the capacitive-coupling scheme in a mi- 
crowave TLR can be used to cool the MR like in the 
case of radiation-pressure cooling in an optical cavity. 



The Hamiltonian of the system reads 



H = hiu'a^a + 



2m 



^^^2 + ;i(£ate— , 



£ ae 



(1) 



The first line describes the free Hamiltonian of the TLR 
and the MR, respectively, with lowering (rising) opera- 
tor of the TLR mode a (a^), and the position (momen- 
tum) operator of the MR x (p) which satisfy [a, a^] = 1 
and [x,p] = ih. The first term in the second line is the 
capacitive coupling between the TLR and the MR. Ac- 
tually, it describes the capacitive energy between them. 
The MR and the TLR are assumed to form a capacitor 
with the capacitance Cg{x) « (1 — x/d) (for small dis- 
placement) depending on the position of the MR along 



III. QUANTUM LANGEVIN EQUATIONS AND 
FINAL MEAN PHONON NUMBER 



The dynamics is also determined by fluctuation- 
dissipation processes that affect both the TLR and the 
mechanical mode. They are taken into account in a fully 
consistent way by the quantum Langevin equations [5l|: 



X — p/m, 



(3a) 



p = -mujlx - 7f,p -I- ^ (2a^a + 1) + ^, (3b) 



d — ~{k + iAo)a + igoax + £ + V'^na 



(3c) 



Here («!„) is the noise operator due to the external 
microwave drive, and ^{t) denotes the quantum Brownian 
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force that the resonator is subject to. They satisfy [5l| 



IV. COOLING OF THE MR 



(aUt)al{t')) = {N + mt~t'), 



(4) 



mm) 



2-K 



dee— (*-*')..(l + coth2^), 

(5) 



where N = l/[exp(?iaja/fcsT') — 1] is the mean number of 
thermal microwave photons of the TLR, fc^ is the Boltz- 
mann constant, T the temperature of the environment, 
and 7ti the damping rate of the MR. For simplicity, we 
have assumed that both the bath correlated to the exter- 
nal microwave drive field and the one connected to the 
MR have the same temperature |43|. We now perform a 
similar calculation as that given in Refs. [3^ [Sa. [stI. |53| . 
The steady-state solution of the quantum Langevin equa- 
tions ([3]) can be obtained by first replacing the operators 
by their average and then setting d(...) jAt — 0. Hence 
we can get the steady-state values as 



(p)=0, (x) 



fipo (|(a>f 



K-l-iA' 



(6) 



where A = Aq — go {x) is the effective detuning. In 
Eq. ©, we can also take |(a)p -I- i ~ since we 

will focus on the case |(a)| ^ 1 which can be achieved 
by controlling the input power of the external microwave 
drive. 

Rewriting each operator as a c-number steady-state 
value plus an additional fluctuation operator, and ne- 
glecting the nonlinear terms (since we have chosen 
I (a) I 3> 1), we obtain a set of linear quantum Langevin 
equations (see Eq. (jAip ) and then solve for the spec- 
trum of the position and momentum of the MR as in 
Refs. Hi, m, In Hi, see Appendix A. 

Using the fluctuation spectra of the MR as given in 
Eqs. (|A10IA11[) . we can define the final mean phonon 
number in the steady state [36| as 



(7) 



2hmLJi, 2h 2 
where the variances of position and momentum are 
1 /■+°° 

{Sr^)^Tr Sr{Lo)dto, {r = x,p). (8) 

This allows us to define the effective temperature T^e as 

(9) 



_ hujh ,1,1 

Kb nl 



In the next section, we will consider the cooling of the 
MR by discussing its final effective mean phonon number 
(or equivalently, its effective temperature) in detail. 



The final effective mean phonon number of the MR can 
be calculated directly by evaluating the integral in Eq. ^ 
numerically and using Eq. ([7]). Alternatively, instead of 
being evaluated directly, Eq. ^ can also be evaluated 
analytically using the approximation scheme described 
in the following. 

The effective mechanical damping rate ■'yf^{uj) = 7f, + 

7ca(w) can be significantly increased^ [jf (w)| ^ jb, 

when |(/o(a)| is very large, see Eq. (|A16p . Let us con- 
sider the most interesting regime when the significantly 
increased effective mechanical damping rate is less than 



the mechanical frequency: |7b*(w)| < oJb, (that is, the 
effective quality factor Q^^ = LUb/\^^^{uj)\ > 1), and 
also less than the decay rate of TLR: |75^(ti^)| < n 
[33I m, [13, [Ell . In this regime, the effective frequency is 
unchanged ujf^{uj) ~ ujb [3a,[5J| according to Eq. (jA15|) . 
and the effective susceptibility is peaked around the 
points iv — ±uj'^^{lu) ~ iwfc. Then one can get an ap- 
proximate expression for the variance 



2m^Lul\jfitOb)\ ' 



(10) 



where the effective thermal noise spectrum 5'(,j(w) and 
the induced noise spectrum 5*^^(0;) are the symmetrized 
parts of Sth{^^) and S'cq(ijj), respectively: 



5*4^ (tj) = /i7bmCi;coth 



2kBT' 



SU^) = i2N + l)hm- 
Similarly, one can obtain 



A^ -I- tj^ 



2A 



■7ca(w)- 



(11) 
(12) 

(13) 



Figure[2]shows the variance of position (fe^) as a func- 
tion of the effective detuning A. The dashed lines corre- 
spond to a numerical evaluation of the integral in Eq. ([5]). 
The solid lines describe the approximate results obtained 
through Eq. (|10p which can be seen to agree perfectly 
with the exact numerical evaluation. We checked that 
this is also the case for the variance of the momentum 

In Eq. ([To]), the induced noise spectrum S'^^{u!b) in- 
creases (heats) the motion of the MR. On the other hand, 
when the effective damping rate is enhanced: |7f;^(wb)| > 
7h, the mechanical motion will reduce, that means cool- 
ing. Mathematically, the cooling effect would dominate 
the heating effect when the effective damping rate is suffi- 
ciently increased. Actually, this is right when the signifi- 
cantly increased effective damping rate is positive for pos- 
itive detuning. However, it is not the case when 7^*^(0;;,) 
is negative and |7jJ*^(wh)| ^ 7b (for negative detuning 
A < 0). That is because the stability conditions, derived 
using Ref. [ssj . are satisfied only for positive detuning 
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FIG. 2: (Color online) Variance of position (Sx^) in units of 
h/mujb as a function of effective detuning A. The dashed 
lines are obtained by numerically evaluating the integral 
in Eq. (ISJ, the solid lines by using the approximate ex- 
pressions Eqs. (|10I13|) . Here, T = 6 x lO^hujb/kB, go = 
3 X lO'^ujt^rnujt/h, cj^ = 2 x W^LJb, e = 2.5 x lO^tJf,, 



76 = 0.25 X 10" 
(lower lines). 



ujb, and K — LUb (upper lines) or «: = 0.2lj6 



[321 152|. I60l|. In fact, a negative effective damping means 
the amplitude motion of the MR will be amplified which 
will lead to an instability [sl, [s^, [111 . 

In what follows, we will focus on the case of positive 
detuning A > 0. According to Eqs. (|7ll0ll3p . one has 



lb + 7ca(t^fc) 



where 



rib 



2hmib^b 



exp{hLUb/kBT) — 1 



(14) 



(15) 



is the initial mean thermal excitation phonon number of 
the MR; 



2N+1 
4cJhA 



1 



(16) 



is the induced mean phonon number due to the capacitive 
coupling between the MR and the TLR. 

As discussed above, the significant reduced value of n[ 
in Eq. (fT4)) can only be obtained when the additional 
damping rate 7ca(wh) (or effective damping rate i^^{ujb)) 
is positive and much larger than the original one (but 
still less than the decay rate of TLR and less than the 
frequency of the MR as discussed before): 



7fc < 7ca(t^fc) < [i^b, k) ■ 



(17) 



This can be satisfied by enhancing the value of \gQ (a)|, 
that is, by controlling the capacitive coupling strength go 
and increasing the external microwave drive power P to 
make |£| large (equivalently, |(a)| will be large). In prac- 
tice, the capacitive coupling strength go would be limited 
by the realistic system, and the external microwave drive 
strength e would also be limited according to the validity 
of the rotating wave approach as we mentioned before. 



FIG. 3: (Color online) Final mean phonon number in the 
steady state nj, vs. effective detuning A and decay rate k of 
the TLR (T = 3 X W'^hLOb/ks, for the other parameters see 
Fig. [2}. 



Here we put the length of the MR I as large as 10 — 100 
fim and the distance between the MR and TLR d as small 
as 1 /xm (see Fig. [T|) in order to get a large capacitance 
Cg which will lead to large go, and fix |e| = Wa/8 for 
all the numerical calculations. Then for a significantly- 
increased effective damping rate, the final mean phonon 
number reduces to 



lb 



Icail^b) 



rib 



(18) 



In order to get the ground state cooling, that is, nl <C 1, 
both rica and ribib/lcai^b) should be much less than 1. 
Especially, if the icai^Jb) is significantly increased enough 
to make 



Ibrib < lcaii^b)nc 
then approaches the limit rica' 



(19) 



(20) 



Now we discuss the possible minimal value of rica by 
discussing all kinds of parameters, e.g., k, A, and N. 
From Eq. (|16p , it is obvious that the optimal value of k 
satisfies the high-quality cavity limit 



and the optimal detuning satisfies A = y'ujf+i- 
Then the corresponding induced mean phonon number 



(21) 



UJb- 



N - 



(22) 



The optimal N needs a sufficiently low initial temper- 
ature of the bath which is limited in practice to the ex- 
perimental dilute refrigerator temperatures. For the su- 
perconducting TLR scheme, its microwave frequency is 
of the order of 2tt x 10^° Hz. For the initial temper- 
ature T > 1 K, ksT > hwa and TV > 1, the ground- 
state cooling of the MR is not possible. Therefore, initial 
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temperature less than 100 mK are required to achieve 
ground-state cooling. 

Our result on the limiting value in Eq. ([22]) is consistent 
with that in other optical schemes except the limit of 
initial temperatures. In the optical cavity case fkUa 2> 
kgT (even at room temperature) and therefore iV ~ 0, 
the optimal value of n^a becomes 



(23) 



which is just the case of resolved sideband coo ling as dis- 
cussed in the optomechanical cooling schemes j22l. [23l [25l. 
[26l . [28j . We would like to point out that these references 
also mention another cooling limit: the Doppler cool- 
ing limit, which is realized in our system when iV ~ 0, 



A = ^/lJf+i^, and > uj^ in Eq. ((T 



2wb 



> 1. 



(24) 



On the other hand, if iV k^/AujI, the induced mean 
phonon number rica in Eq. (|22p becomes rica — *■ N. In 
the classical limit when the initial temperature is so high 
that N « ksT/^hLUa) > 1, one has 



7c ft 



rib 



rib 



UJb ' 



(25) 



which is also given in Ref. [3l| . The Doppler cooling limit 
in Eq. ([M]) and the classical cooling limit in Eq. (P5|) pre- 
clude ground state cooling. We will focus on the resolved 
sideband cooling in this paper. 

The final mean phonon number is plotted as a func- 
tion of the effective detuning A and the decay rate k of 
the TLR in Fig. [31 It is clear that one can obtain a 
significant suppression of the mechanical motion of the 
MR in the positive detuning range A ~ uJb- The opti- 
mal cooling is obtained for k? <^ lo^, which agrees with 
both the above analysis and that in other treatments of 
radiation-pressure cooling [H, [H, [H, [H, \^ . 

Physically, as discussed in the back-action optome- 
chanical cooling schemes in optical cavities (Refs. [Sl, 
[3^ . [3^ . [37I Hoj), the external driving microwave is scat- 
tered by the "TLR -I- MR" system mostly to the first 
Stokes sideband [ud — Wb) and the first anti-Stokes side- 
band (iOd + Wb). The generation of an anti-Stokes photon 
will cool the MR by taking away a phonon of the MR. On 
the contrary, the generation of a Stokes photon will heat 
the MR by creating a phonon. When the effective de- 
tuning A > 0, the microwave field of the TLR (with the 
frequency Wq = w^-f Aq ujd + A) interacts with the first 
anti-Stokes sideband (cj^-l-a^h) more than it interacts with 
the first Stokes sideband {ojd—i^b), and cooling will occur. 
This is the physical reason why the positive effective de- 
tuning (A > 0) will lead to cooling. In the high-quality 
cavity limit k < Ub, the anti-Stokes (Stokes) sideband 
is resolved, and the corresponding cooling (heating) pro- 
cess is prominent. Especially, for the optimal effective 



detuning A « w^, the frequency of the TLR is resonant 
with that of the anti-Stokes sideband, which will appar- 
ently lead to optimal cooling. This physical discussion is 
consistent with the calculation presented above. 

Figure [3] suggests there is a finite optimal value of 
K for a fixed effective detuning. That is because one 
should have both a small value of iica and a large effec- 
tive damping rate 'jcai^b) in order to get strong cooling: 
K should not be too large, since Uca depends somewhat 
on the value of k^/w^; k should not be too small, since 
7ca(wb) — > when k — > 59]. In the coohng process, the 
thermal energy of the MR is mainly first transferred to 
the TLR, and then leaks out of the TLR through the bath 
coupled to the TLR. When the decay rate of the TLR is 
too small: k ^ 0, the energy leakage out of the TLR is 
too weak, and one could not obtain a strong cooling. 



0.002 




FIG. 4: (Color online) (a) Optimal final mean phonon 
number obtained by numerically evaluating the integral in 
Eq. ([51 (dashed line) or using the approximate expression in 
Eqs. (|10ll3p (solid line) as a function of k at the optimal ef- 
fective detuning A — LJb- (b) Logarithm of the ratio of the 
corresponding effective damping rate ^^^{LJt) to k as a func- 
tion of K. Here, T — 3x IG^hujb/kB- For the other parameters 
see Fig. [1 

We would like to emphasize that the results shown 
in Fig. [3] are based on the approximate expressions 
Eqs. (|10ll3p . where the condition of the so-called weak- 
coupling limit [33, SI] has been assumed, that is, the 
effective damping rate of the MR should be always less 
than the decay rate of the cavity and less than the fre- 
quency of the MR l7fc*^(ti'h)l < k, Wft. Normally the weak- 
coupling is satisfied but not in some special cases. In 
Fig. m^b) , the weak-coupling condition is violated when 
n/ujb < 0.1 at the optimal effective detuning A = cub- 
Beyond the weak coupling limit, Fig.UI^a) shows that the 
approximate treatment through Eqs. (|10ll3p ceases to be 
valid. Then one should discuss the cooling, e.g., effective 
mean phonon number in Eq. ([7]), by using the numerical 
evaluation of the integral in Eq. ([5]). But going beyond 
the weak-coupling limit, the contribution from the posi- 



6 



tion variance is not equivalent to that from the momen- 
tum variance any more [36| . In other words, the energy 
equipartition is not satisfied. That means it is hard to 
define an effective temperature since it is not in a strict 
thermal state. 

According to the above analysis, both the high-quality 
cavity and weak-coupling limit should be satisfied, so 
the optimal decay rate of the TLR is better taken to 
be K w O.lwb for the typical parameters in Fig. [4] The 
weak-coupling condition depends only weakly on the ini- 
tial temperature T and the original damping rate of the 
MR 7b in the cooling process. In what follows, we will 
consider the optimal decay rate at k « O.lwf, for different 
parameters T and 7b, for which the weak coupling limit 
is always satisfied. 

In Fig. [HI the ratio of final effective temperature Toff 
to bath temperature T is plotted as a function of the 
effective detuning A for the optimal k « O.lwfc. Appar- 
ently, here the weak coupling limit is satisfied (according 
to the above analysis in Fig. [3]). For initial tempera- 
tures T = 100, 30, and 10 mK, the corresponding initial 
mean phonon numbers are — l/lexp^hujb/kBT) — 1] 
~ kBT/hujt ~ 3300, 980, and 330, with the final mean 
phonon number nl « 1.6, 0.5, and 0.16, respectively. It 
is obvious that a significant cooling of the MR is ob- 
tained and lower initial temperatures will generally lead 
to better cooling. For an initial temperature T = 10 mK, 
which can be realized experimentally by using a dilution 
refrigerator, the MR (with the frequency Wfc ~ 4 MHz) 
can be cooled close to the ground state since the final 
mean phonon number w 0.16 < 1. 




0.6 0.8 1 1.2 1.4 

A/co^ 

FIG. 5: (Color online) The final mean phonon number vs. 
effective detuning A for three initial temperatures: T = 10 
mK (dotted lines), T = 30 mK (solid lines), T = 100 mK (dot- 
dashed lines). Here, k — O.lojj,, m = 1.5 x 10^^^ kg, = 4 
MHz, 76 = 0.25xlO"*tJ6 (equivalent^, Qt = Ldb/jb = 4x10"). 
For the other parameters see Fig. [5] . 

For an initial temperature of T = 10 mK as in Fig. [5l 
the final mean phonon number would in principle be 
(in Eq. p^ ). which is much less than that obtained in 
Fig. El Uca = N + k'^/AujI « 0.0025 < « 0.16. This is 
because n\ rica only when the condition in Eq. is 
satisfied. Unfortunately, it is not the case for the parame- 
ters in Fig. O A possible way to approach this condition 



is to increase the quality factor of the MR. In Fig. [6l 
the final effective mean phonon number is plotted as a 
function of the effective detuning A for different quality 
factors of MR Qh (= tOb/nb): Qb = 4 x 10'' (typically, 
see Ref. [H]); Qb = 10^ 4 x 10^ 10^ (expected in the 
near future). The corresponding minimal w 0.16, 0.06, 
0.02, 0.01. One can find the cooling is better for a higher 
quality factor of the MR. 
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FIG. 6: (Color online) The final effective mean phonon num- 
ber n[ is plotted as a function of A for different quality factors 
of the MR: Qb = 4 x lO", 10^ 4 x 10^ lO" (from up to down), 
the corresponding damping rate 75 = 100, 40, 10, 4 Hz. Here, 
T = 10 mK. For the other parameters see Fig. [S] 

The cooling discussed above can be measured by a ho- 
modyne detection method like that given in the scheme 
of cavity-assisted radiation-pressure cooling of a MR 
[3^ . H, 113, mi . The motion of the MR can be detected 
by monitoring the output microwave signal (e.g., the field 
phase quadrature) of the TLR (as seen in Fig. [1]) since 
the measurement of the output spectrum corresponds to 
a faithful measurement of the MR motion [47] . 

V. CONCLUSION 

We have found that a MR with frequency ojt ^ 27rx 10^ 
Hz can be cooled close to its ground state when it is cou- 
pled to a typical TLR (wa 27r x 10^° Hz). Actually, by 
considering the optimal parameters in this scheme, that 
is, assuming the high-quality cavity limit (f? ^ 
a positive optimal effective detuning (A k, cjb), a low 
initial temperature (e.g., T = 10 mK in order that 
N K, 10"^"^ ~ 0), a high quality factor of the MR 
(Qb = '-^b/lb ^ 10^), and both strong external input 
microwave drive power P and strong capacitive coupling 
strength to get the significantly increased positive ef- 
fective damping rate (7^ ^ ^l^{ijJb) ~ -fca{i^b)), we find 
that resolved sideband cooling of the MR occurs. The 
possible minimal value of the final effective phonon num- 
ber could approach the induced mean phonon number: 
n-l Uca ~ >^ l^'^l 1- Moreover, one should also 
consider the condition of weak-coupling limit, that is, 
the significantly increased effective damping rate 7^*^(0;^) 
should be less than both cub and k (in the high-quality 
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cavity limit, one only needs "fl^{ujb) < k). This condi- 
tion requires that k must not be too small though lower k 
will lead to lower rica- As shown in Fig. [2] and its discus- 
sion, there will be an optimal range of k. For the typical 
parameters in this cooling scheme, we take k ~ O.lcJb 
(though this is not the optimal result in general). We 
find that a MR with Wf, = 4 MHz can be cooled close 
to its ground state with the final effective mean phonon 
number in the steady state: nl « 0.f6 (for a typical qual- 
ity factor = 4 X f 0^) or nl ~ 0.01 (for a high quahty 
factor Qh ~ 10®) by using the resolved sideband cooling 
scheme when it is coupled to a driven TLR (with the 
frequency = 8 x 10^" Hz). 

We would like to stress the condition in Eq. p9| , which 
can lead to — > Uca- As pointed out in the discussion 
of Figs. ([SI [6|), this is not always satisfied. A possible 
way to approach this condition is to increase the quality 
factor of the MR. For example, if the quality factor of 
the MR is high enough (e.g., Qb > 10^), one would have 
the optimal — > n^a ~ 

/AujI = 0.0025, for which the 
MR is cooled much closer to the ground state. 

The back-action self-cooling scheme presented here 
is similar to the optical-cavity-assisted cooling scheme 
[H, [IJ- In both cases, the MR can be cooled close to 
its ground state using resolved sideband cooling which is 
possible in the limit of a high-quality cavity. But it seems 
that this limit (e.g., k = O.luJb, for = 4 x 10® Hz) is 
easier to reach in the microwave TLR than that in the op- 
tical cavity. In the case of a TLR (typically = 8 x 10^° 
Hz), quality factors of Qa = oJa/n = 2 x 10^ have been 
seen in experiments [13, [fi^l- However, in the case of an 
optical cavity {ua ~ 10^^ Hz), the corresponding qual- 
ity factor should be Qa = LJa/K ^ 2.5 x 10^, which is 
hard to achieve since the typical cavity quality factor is 

To conclude, we have studied the self-cooling of a 
mechanical resonator that is capacitively coupled to a 
transmission-line resonator. The discussion was based 
on a linearized quantum Langevin equation. The cool- 
ing method presented here is similar to the self-cooling 
of a MR coupled to an optical cavity by radiation pres- 
sure. By using the optimal parameters discussed above, 
the MR can be cooled close to its ground state in the 
high-quality cavity and weak-coupling limit. 
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APPENDIX A: EQUIVALENCE TO 
TIME-DEPENDENT SECOND-ORDER 
PERTURBATION THEORY 

The linearized quantum Langevin equations read 

6x = 6p/m, (Ala) 
5p = —niLulSx — jbSp + hgQ{6a) (a) + h.c.) + ^, (Alb) 
5a — — (k + iA)(5a -f ig^ (a) Sx + \/2Kai„. (Ale) 

To solve these equations, we define the Fourier trans- 
form for an operator u (u — Sa, Sx, 5p, ai„, ^) 



u{t) 



1 



/2^ 



e*'^*u(a;)dtj. 



and for its Hermitian conjugate (if any) 

\ /• + 00 



u\t) 
which lead to 



(A2) 



(A3) 



(A4) 



^{9)i{uj) ) = h-ibmuj{l + coth ^rr^)5{n + uj). (A5) 
/ 2kbJ 

After solving the linear quantum Langevin equations 
in the frequency domain, we obtain 



6x{uj) 



B{u:) 



5p{uj) — iu!m6x{Lu), 



(A6a) 
(A6b) 



where B{uj) — m{uj'^ ~ uj"^ + ijb^) [(« + i^)^ + A^] — 
2h\go{a)\'^ A, and C{ui) = hV2K,go{a)[K + i{uj + A)]. 

To calculate the effective temperature of the MR, we 
define the fluctuation spectra of position and momentum 
of the MR, which are given by the following 
correlation function: 



Sp{uj) = 



-foo 



e-"^^ {Sx{t + T)Sx{t))^dT, (A7) 
e-'^^ {Sp{t + T)dp{t))^dT. (A8) 



Here, {■■■)s denotes the steady-state average. Equiva- 
lently, Sx,p{(^) can also be defined as 

{SP{n)SP{Lo))^ := Sr{Lo)S{n + uj), {r = x,p). (A9) 



According to Eqs. (|A6p . the spectra of the MR can be 
written as 

Sxioj) = |Xcff(^)|' [Sthiio) + S,ai^)] , (AlO) 



Spiuj) = (wm) Sxiuj), 



(All) 
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where 



Sthi'-^) — hjhmujll + coth{hLu/2kBT)] 



(A12) 



is the thermal noise spectrum due to the Brownian mo- 
tion of the MR; and 



Sca{uj) 



{N + l)\C{u;)\-' + N\C{-aj)\' 



|(k + zcj)2 + A2|2 



2ujA 



K + + A2|2 



(A13) 



is the induced noise spectrum due to the capacitive cou- 
pHng to the driven TLR The effective susceptibiUty is 
defined as Xoff(w) = [(k + ioj)"^ + A^] /B{ijj) and can be 
simplified to 



1 



XcSil^) = r ^ 

where the effective frequency of the MR is 



(A14) 



^^,_ 2n|go(«)rA(.2,.. + A2) ^^^^^ 



(k2 -tj2 + A2)7,„(tj) 



2k 



and the effective damping rate ^f^i^^^ = Ih + lcaij^') with 
the additional term 



4;i|.go(a)| kA 
TO |(k + iujY + A2|^ 



(A16) 



resulting from the capacitive coupling. 

According to the definition of the additional damping 
rate in Eq. (|A16p . the induced noise spectrum 5'ca(w) in 
Eq. (jAfSj can also expressed as 



(2iV+l). 



2A 



- LO 



(A17) 
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